IGC-09/5-3 
PI-QG-134 

A geometric perspective on singularity resolution and uniqueness 

in loop quantum cosmology 

Alejandro Corichi^'^'B and Parampreet Singh^'Q 

^Instituto de Matemdticas, Unidad Morelia, 
Universidad Nacional Autonoma de Mexico, UNAM-Campus Morelia, 
A. Postal 61-3, Morelia, Michoacdn 58090, Mexico 
^Center for Fundamental Theory, Institute for Gravitation and the Cosmos, 
Pennsylvania State University, University Park PA 16802, USA 

^Perimeter Institute for Theoretical Physics, 
31 Caroline Street North, Waterloo, Ontario N2L 2Y5, Canada 

Abstract 

We re-examine the issue of singularity resolution in homogeneous loop quantum cosmology from 
the perspective of geometrical entities such as expansion rate and the shear scalar. These quantities 
are very reliable measures of the properties of spacetime and can be defined not only at the classical 
and effective level, but also at an operator level in the quantum theory. Prom their behavior in 
the effective constraint surface and in the effective loop quantum spacetime, we show that one can 
severely restrict the ambiguities in regularization of the quantum constraint and rule out unphysical 
choices. We analyze this in the flat isotropic model and the Bianchi-I spacetimes. In the former case 
we show that the expansion rate is absolutely bounded only for the so called improved quantization, 
a result which synergizes with uniqueness of this quantization as proved earlier. Surprisingly, for 
the Bianchi-I spacetime, we show that out of the available choices, the expansion rate and shear 
are bounded for only one regularization of the quantum constraint. It turns out that only for this 
choice, the theory exhibits quantum gravity corrections at a unique scale, and is physically viable. 
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I. INTRODUCTION 



One of the main challenges that faces a potential quantum theory of gravity is to resolve 
the classical singularities that are present in general relativity (GR), and in particular the 
initial singularity, or the big bang. Loop quantum gravity is a candidate for such a 
theory. In particular, when loop quantization methods are applied to isotropic homogeneous 
cosmologies, in what is known as loop quantum cosmology (LQC) [?] the resulting theory has 
shown interesting promise in this enterprise. For, when models with simple matter content 
are considered, LQC yields a well defined quantum theory for which the big bang/big crunch 
gets replaced by a quantum bounce. These results first derived for a spatially flat FRW 
universe with a massless scalar field [1, 0, [1] were generalized to include the k = ±1 cases 
j^, 0], a cosmological constant ^ and also a 0^ inflationary potential [SS]. Within the k = 
Robertson- Walker model and massless scalar field, the model could be solved exactly by 
choosing volume as the lapse [loj. This enabled us to prove robustness of various results 
which were found previously through numerical simulations. These include: (i) quantum 
bounce for a dense subspace in the physical Hilbert space and (ii) a supremum on the 
expectation values of the energy density operator which coincides with the critical energy 
density at which the bounce occurs. Further, using this model it is possible to gain insights 
on the nature of the universe before the quantum bounce. It turns out that a universe 
which is peaked on a classical trajectory at late times after the bounce is constrained by the 
properties of the quantum constraint and observables to be peaked on a classical trajectory 



before the bounce at early times [llj. For an up to date review of these results see [l2j. 

It is well known that a process of quantization may result in various quantization ambigu- 
ities. More importantly, regularization of the constraints can be performed in inequivalent 
ways of which not all may lead to a physically consistent and viable theory. One thus re- 
quires a criterion to eliminate the choices which are unphysical. For the k = isotropic 
model these criteria were put forward in Ref. where different loop quantizations were 
compared and put to test by requiring a set of reasonable consistency conditions. It was 
shown that there is only one possible regularization of the constraint operator for which the 
following conditions are satisfied [l3| : (i) independence from the freedom underlying fiducial 
structures, (ii) a well defined critical energy density at which bounce occurs and (iii) that 
at small spacetime curvatures the theory should approximate GR. The only loop quanti- 
zation in the flat isotropic case which satisfies these requirement is the so called 'improved 
quantization' 

In these studies the use of effective dynamical equations turns out to be very helpful. 
Effective equations are ordinary Hamilton equations of motion found by means of an effective 
Hamiltonian Ties- This Hamiltonian (that depends explicitly on h so it can not be seen as 
purely 'classical') has been found using geometrical methods in quantum mechanics [13] and 
has been shown to be responsible for the dominant term to the canonical path integral of 
the theory [l^. In a recent systematic study of the resulting 'effective spacetimes', it was 
shown that, for the k = case, they are always non-singular and there can exist no strong 
curvature singularities [l6|. Thus, in a precise sense and for k = isotropic case, loop 
quantum cosmology yields a generic singularity resolution. 

Even when these results are rather encouraging towards the ultimate goal of proving 
that quantum effects generically resolve the classical singularity, one still needs to overcome 
several challenges. For instance, when considering anisotropic models, like a Kantowski- 
Sachs and Bianchi I cosmologies, the notions that were rather useful for isotropic models 
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already show some limitations. When we depart from the isotropic models, the first difficulty 
is that the curvature tensor not only has a Ricci component but also its Weyl part is non- 
zero. This means that the classical spacetimes can develop 'Weyl singularities', in which the 
curvature invariants build from the Weyl part of the curvature blow up (while allowing for 
the Ricci part to remain well behaved). Thus even if quantum effects bind the Ricci scalar, 
the spacetime may still contain physical singularities.^ These considerations bring to the 
forefront the need for an unifying principle that allows us to approach the issue of generic 
singularity resolution in the quantum theory. In other words, even at this level one is pressed 
with the question of what criteria should one use to say that singularities are resolved by 
the quantum theory. 

The question in itself poses its own problems. Singularities in the classical theory are 
understood in terms of concepts such as the inextendibility of geodesies on the spacetime [13] ■ 
However, the mere notion of geodesies assumes the existence of a sufficiently differentiable 
metric for all points on the manifold, including those that are in the region where 'quantum 
effects' should dominate. Quantum gravity is expected to introduce some 'fuzziness' near 
the singularity, which itself means the notion of spacetime, and geodesies, should stop to be 
valid and be replaced by some 'quantum foam'. How can we distinguish then the situation 
where singularities are resolved from those where they are not, if the quantum state is 
not 'semiclassical' in the strong gravity regime? In other words, is the singularity resolved 
just because the state is not semiclassical near the would-be singularity? Several proposals 
have been put forward in this direction. These include the idea of singularity resolution 



via quantum foam effects |18l . |19| . For a recent review of different such scenarios, from 
the perspective of Wheeler-DeWitt theory, see for instance |2Ci]. From our viewpoint this 
question can only be meaningfully addressed via physical observables. This implies, in the 
usual language of constrained systems, analyzing the properties of Dirac observables. It 



perspective that singularity resolution was established for the FRW models in 



is from this 

LQC [3, @, [7, 10, 21], when the relational observable defined by the matter energy density 
'at internal time T' was shown to be absolutely bounded. 

The purpose of this paper is to revisit the problem of singularity resolution for homoge- 
neous cosmologies from a geometric perspective. Here we put forward the viewpoint that 
the relevant quantities to consider are invariant geometrical quantities that are in the clas- 
sical setting responsible to define when singularities arise. In particular, we shall pay close 
attention to well motivated scalars associated to the evolution of cosmological observers such 
as the expansion parameter and the shear^. As we shall argue, these geometrical quantities 
not only are useful for considering the evolution as given by the effective Hamiltonian, but 
can indeed be extended to the quantum realm where the classical interpretation might not 
be valid. This allows for a quantum version of singularity resolution. 

Another aspect that we shall consider here concerns the issue of choosing a unique quan- 
tization for the constraint operator of the theory. In the case of isotropic cosmologies, math- 



ematical and physical consistency condition uniquely selects a quantization [13|]. A natural 



^ The situation in the anisotropic models is thus more complex than the isotropic model where a bound 
on the energy density of the matter field translated to the bound on spacetime curvature for all matter 
(except whose pressure diverges at finite density) and the only singularities which may exist are harmless 
weak curvature type In the anisotropic model, the bound on energy density thus captures little about 
the fate of singularity. 



As is well known, expansion and shear play a prominent role in the classical singularity theorems 17 1. 
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question is whether this quantization also gets selected from the geometric perspective of 
singularity resolution. As we shall see, in the case of isotropic cosmologies, the answer is in 
the affirmative. The quantization which is naturally selected is the one resulting from the 
'improved quantization' of LQC [1, With this in mind, we analyze two quantizations of 
the Bianchi I model proposed in literature, for which one expects to recover the isotropic sec- 
tor. Both these quantizations are motivated by a generalization to anisotropy of the unique 
'71' quantization of the isotropic sector, but are inequivalent. As we show here, requiring that 
the theory has a well defined Planck scale, and therefore resolves the singularity from our 
geometric perspective, picks out one of these quantizations. This is the quantization which 
has been recently developed in detail by Ashtekar and Wilson- Ewing [22|. Our analysis thus 
leads to a lesson for the anisotropic model, which is: even if inequivalent quantizations might 
be expected to lead to the correct isotropic limit, a detailed analysis shows that this might 
not be so. Therefore, not all of them are viable. In fact, for the Bianchi-I cosmologies our 
results point towards a similar sort of uniqueness as for the isotropic LQC. 

The structure of the paper is as follows. In Sec. [Ill we give some preliminaries regarding 
the behavior of cosmological observers within homogeneous cosmologies, including the mod- 
ifications brought in by the effective dynamics of LQC. In Sec. IIIII we review the isotropic 
FRW model from this perspective and show that the 'improved' quantization in tailored to 
bound the expansion of the geodesies of cosmological observers and therefore gets selected 
as a viable choice. Section IIVI is devoted to the Bianchi I model for which we compare two 
existing quantizations in view of our geometrical considerations. We end with a discussion 
and an outlook in Sec. El 



II. HOMOGENEOUS COSMOLOGIES AND OBSERVERS 

With the purpose of having a self-contained exposition, in this section we recall some 
basic notions needed to treat cosmological observers in homogeneous settings. We include 
a discussion of the corresponding modifications that emerge from the effective dynamics of 
loop quantum cosmology for the isotropic case. 

Let us first recall the main features of the congruence of cosmological observers and the 
way their expansion is related to other geometric invariants. If denotes the congruence 
of cosmological observers, then the quantity 

Bab:=S/b^a, (1) 

is symmetric and corresponds to the Lie derivative of the metric in the normal direction 
(and for = 1 this is just the proper time derivative), and also to the extrinsic curvature 
Kab of the slicing by homogeneous spaces. The trace of this quantity is 

e:=g'^^B,, = K . (2) 

That is, the expansion of the congruence of cosmological observers corresponds to the trace 
of the extrinsic curvature. Given that the function 9 contains the information about the 
expansion of the cosmological observers, it is a natural measure of when singularities form, 
namely when the expansion blows up. The evolution equation for the expansion rate 9 
relates it to the shear scalar cx^ and the spacetime curvature: 

^ = -\e'-a^''aab-Rabee (3) 
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where the contribution from the term corresponding to vorticity vanishes due to the existence 
of homogeneous shces. This equations holds in general for homogeneous cosmologies for 
which a preferred set of observers form the congruence 

For classical GR, if matter satisfies the strong energy condition then Rab^"' > 0, and 
hence the expansion rate of an initially converging congruence diverges in a finite proper 
time. The classical constraint (density) relates the extrinsic curvature Kab, the spatial Ricci 
scalar (which is zero in our case) and the energy density (p), 

K''''Kab-K'^ = lQ7rGp (4) 

which can be rewritten in terms of 6 and cr^ as 

16nG p = 9^ - (t\ (5) 

For the isotropic spacetime, the divergence in expansion rate or the caustic is equivalent to 
the divergence of energy density and also that of curvature invariants such as Ricci scalar: 

R = SttG (1 - 3w)p . (6) 

Here w is the equation of state of the matter source: w = P/p, where P denotes the 
pressure of the matter content. For the anisotropic spacetime, a caustic may be formed 
from divergence either of energy density and/or the shear scalar. In the former case the 
curvature singularity is Ricci curvature type and in the latter it is Weyl curvature type. 

If we now consider the Hamiltonian theory, Eq. can be rewritten in terms of an 
arbitrary choice of canonical variables, denoted by (p, q) . Let us represent these function by 
6 = 6{p, q) and {p, q) . It is important to stress that these functional relations are 

already using some dynamical information (relating for instance, extrinsic curvature with 
time derivative of the three- metric, which are part of Einstein's equations in canonical form). 

In the effective dynamics of LQC though Eq.Q holds, the relationship between expansion 
6, shear and p is modified because the classical constraint (jlj) is no longer valid. In loop 
quantization, the classical constraint in terms of Ashtekar variables i.e. the connection A 
and triad E is first expressed in terms of holonomies of connection and flux of triads and then 
quantized. For the fiat isotropic model and the Bianchi-I spacetime, the only modifications 
to the classical constraint arise in the gravitational part.'^ It implies that eq.([5]) the scalar 
constraint (and therefore, Einstein's equations) is modified as 

lQnGp = fie,a), (7) 

where f{9, a) is a function determined by the quantization prescription. In the effective 
dynamics, in order to recover general relativity, one expects that f{6,a) h-> | 0^ — cr^ in 

the appropriate regime (when O'^^a'^ <^ V^pi)- ^^i^ ^^^e then we expect that the 

dynamical equations associated to the effective Hamiltonian will approximate Einstein's 
equations also in this regime. 



^ For models with non-vanishing intrinsic curvature, modifications to the inverse scale factor operator affect 
matter energy density. If the intrinsic curvature vanishes, such modifications have little physical meaning 

a. 
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In order to illustrate the way key equations get modified, let us consider the case of a 
k = Robert son- Walker metric where all expressions can be explicitly computed (recall that 
in this case 5"^ = 0). In this case, the loop quantum effective Hamiltonian for the 'improved 
quantization', when described in canonical variables {f5, V), is given by 



3 sm'jXP) 



where V denotes the physical volume and A = A^/Sn^yip^ with 7 = 0.2375 22|. The value of 
A is determined by the procedure of regularization of the Hamiltonian constraint in LQC 
Thus, for LQC, the function f{6, a) in the flat isotropic model is 

a sin^(A/3) , . 

where [3 is canonically conjugate to the volume V and satisfies V} = inG'j. (On the 
classical solutions f3 is related to the scale factor as f3 = 'ya/a). The quantum gravity 
parameter A has dimensions of length and is of the order of Planck length. This parameter 
fixes the scale near which quantum gravity become significant. In fact, when A/? -C 1, we 
recover the classical constraint 9 = /?^/7^ = ^y^p. 

In order to evaluate the expansion factor in the effective LQC we first use Hamilton's 
equation: V = {V,H} and then compute 9 := V/V. For the effective Hamiltonian Hes, we 
get, 

9 := ^ {V,n,s} = ^ smiX/3) cos(A/?) = sin(2 X/3) (10) 
V A7 2A7 

Note that we have retained the use of 9 to denote the expansion as seen by the cosmological 
observers, in the underlying spacetime. In the classical limit, A ^— 0, we recover the relation 
^ ~ 3/5/7 = ^.^ Using the effective constraint ([8]) along with above Hamilton's equation we 

9 3 V Peri./' ^ ' 

with pcrit = 0.41ppianck the Critical density for which the Hubble parameter vanishes H = 
and a bounce occurs. At classical scales i.e. p ^ pcrit, 9^9 and the above equation 
approximates eq.([5]) for the isotropic model. 

Using the Hamilton's equation for (3, it is straightforward to evaluate 9 + 9"^/ 3 = 3d/a 
and find, from Eq. Q 

Rah = ^TrGp f 1 - 4 + UtiGP (i-2 . (12) 

V Pcrit J \ Pcrit/ 



Note also that even when 9 « 3/3/7 when A/3 <C 1, it is different from the expression 9' := 3sin(A/3)/7A 
as it appears in the constraint ([S]). This difference between the quantity that replaces 9 in the modified 
constraint ([5]), and the expansion as computed using the dynamical equations, is what allows to have an 
'effective Friedman equation'. This is qualitatively different to the classical equation (equivalent to the 
classical constraint (O), when written in terms of density, by using the identity sin^(2A/3)/4 = sin^(A/3)(l- 
sin2(A/3)). 
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This implies that in effective dynamics of LQC, Rah^""^ can become negative near the Planck 
scale, even for matter which satisfies the strong energy condition. It is also illustrative to 
rewrite Raychauhuri equation as, 



As an example, for the massless scalar field, where (p = P), the quantity Rabi"'^^ becomes 
negative when p > 2/5 pcrit- Since Rabi"'^^ > is violated in this regime, there is a repulsive 
force that out- weights the gravitational attraction, in such a way that the expansion reaches 
its maximum absolute value when p = 1/2 pent and then \6\ begins to decrease until it 
reaches zero at the bounce, when p = pcrit- Note also that this qualitative behavior of the 
expansion is valid when the matter satisfies the null energy condition. Therefore, caustics 
do not form in the effective spacetime of fiat and isotropic LQC. In view of this, it is not 
surprising that this spacetime is devoid of any physical singularities |il6i] . 

III. ISOTROPIC /c = FRW 

The flat, isotropic and homogeneous Robertson- Walker metric is one of the simplest cos- 
mological model which can be quantized in LQC. Even though simple, this model offers 
useful insights on the regularization of the Hamiltonian constraint and restrictions on quan- 
tum ambiguities for models with more degrees of freedom and even the full theory. It is 
well known that out of the available "loop quantizations" for this model, only the so called 
"improved quantization" [3] leads to a physically viable description. Unlike other regulariza- 
tions of the Hamiltonian constraint, the one described by the (/?, V) variables (as in Sec. [Tll) 
leads to a well defined Planck scale and approximates GR at low spacetime curvature for 
matter satisfying the null energy condition jl3|. 

In the following we re-analyze the quantization of the fiat FRW model for different 'loop 
quantizations' and consider the behavior of the corresponding operator for the expansion 
rate of cosmological observers in all these theories. Perhaps not surprisingly we find that, 
even when many such quantizations might display a 'bounce', there is a unique choice for 
which the corresponding operator is bounded, and therefore, when interpreted in terms of 
the effective description, the resulting bounce occurs at a unique 'Planck scale'. This choice 
corresponds to the improved quantization of j3] ■ In all other choices, including the so called 
'po quantization' (and all lattice refining models [23|), the corresponding operator is not 
bounded. As expected, the quantum gravity scale in this model turns out to be arbitrary 
and the theory lacks any predictive power. In this case, when the effective description 
is studied, the maximum value of the expansion (or acceleration) a geodesic congruence 
encounters can be arbitrarily low, well below the 'Planck scale' 9 k, l//pi. One is then led 
to conclude that there is a unique quantization prescription for which singularity resolution 
is well defined on the full space of states of the theory. We find therefore an interesting 
interplay with the results of [3]. 

Let us recall the most basic facts about the k = models. In the ordinary geometro- 
dynamical description, the geometric configuration variable is the scale factor a and its 
corresponding conjugate momenta Pa = — 3Vo a a/(47rG), with Vq the fiducial volume of the 
region under consideration (for simplicity we shall consider a closed topology for the 




(13) 
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spatial manifold).^ When using the basics variables that are employed in loop quantum 
gravity, namely, densitized triads and connections, it is convenient to consider the variables 
\p\ = Vo^^ a? and c = 7 Vo^^ a, valid on solutions, for which we have canonically conjugate 
variables (c, p), satisfying 

r T 87rG'7 , . 

{c,p} = ^— . (14) 

Starting from these variables, one can consider a rather general possibility of variables, 

(l-m) 

Ps = cp"', 9 = ^ (15) 

1 — m 

obtained from (c, p) by a canonical transformation, for which the effective Hamiltonian 
constraint becomes |l3| 

3 ^|^;^^(^^^ X x(l-4m)/(2(l-m)) _ ^ (-.^s 

^^53^ ((1 m)g) -H^^,, (16) 

From here one can compute the expansion rate 6 = V/V, and get 

3 g 



2{l-m)g 

9(1— m)~^'^ (2m-l) 

= sin(Ap,P,)cos(Ap,P,)^7"^(--) . 

It is now straightforward to see that the absolute value of the expansion will reach a max- 
imum value |6'|max before decreasing and becoming zero (when cos{Xp^Pg) vanishes and the 
bounce occurs). It is also easy to see that the maximum value, on a given trajectory, also 
depends on the phase space variables. In particular, the value of 16*1 max is determined by the 
initial conditions for the matter content except for one case. This is for the choice m = —1/2, 
which yields Pg = [3 and g = 2V/3, and therefore corresponds to the 'improved JI quantiza- 



tion' [4, Il3[. Only in this case \9\ is absolutely bounded, on the constrained surface, by a 



sup- 



phase space independent number | 6 

In order to appreciate the significance of this result, let us consider a particular example, 
namely the so called '/xq' quantization with a massless scalar field (for a discussion of other 
inadequacies of this quantization see j3, [IS])- In this case m = 0, and |^|max is finite for 
any trajectory, but it depends on P^, a Dirac observable that labels the trajectory on phase 
space. Since 9 is bounded on each trajectory (though by dependence on P^), one may 
naively conclude that the theory achieves singularity resolution. As we shall now argue this 
conclusion is incorrect. This can be seen in two steps. First, if we only consider the effective 
equations, and pretended that there is an underlying spacetime at all points of the 'evolution' 
on phase space, the scale (with dimensions L~^) at which quantum gravity effects become 
important can have any value, since the scale depends on P^ which is completely free. This 
means that quantum effects can kick in at an arbitrarily small scales in contradiction to 



^ The results which foUow are not affected by any other choice including that of non-compact R"^ manifold. 
In the latter case the following discussion is supplemented with subtleties of the freedom of the rescaling 
of the fiducial cell which is discussed in detail else where . 
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various experiments, or at an arbitrary scale well beyond the Planck scale. Clearly this 
scenario is unphysical. 

But, is this behavior an artifact of considering the effective equations to be valid at all 
regimes? Can we forgo this pathological possibilities when we consider the full quantum 
theory without going to the effective description? In other words, can we still have quantum 
singularity resolution in the '/xq scheme'? Let us now see that, when the quantum theory 
is considered and the question is carefully analyzed, the answer is still in the negative. For 
that, let us consider in the quantum theory, for arbitrary m, an operator | 9\ associated to 
the expansion, where we 'promote' the (effective) classical expression to an operator in the 
underlying geometrical sector of the Hilbert space of the theory. Clearly, this quantity has 
the standard geometrical interpretation for those states that are semiclassical and whose 
dynamics approximate that of a geometry satisfying classical equations. An important 
feature is that, being a well defined operator, it can also be defined on those state with no 
semiclassical interpretation. If we now compare the corresponding operators for all possible 
choices of m, we are led to conclude that it is only for the improved quantization that the 
operator is bounded in the physical Hilbert space of the theory. Let us be more precise. 
The function 9 that represents expansion for this particular quantization, is bounded on the 
constraint surface. Therefore, the corresponding operator on the quantum theory will have 
a bounded spectrum. This is precisely realized in the exact quantization of [lo|. This result 
is, in this simple isotropic case, rather similar to the behavior of energy density that was 
found to be bounded in this same quantization 



These considerations suggest the following conclusions. First, one could have as a nec- 
essary condition for singularity resolution in the full quantum theory that the operator | 9\ 
be a bounded operator. In the = case, this immediately implies that this condition will 
only be satisfied by the improved quantization, pointing again to a uniqueness result along 



the lines of [13|. Secondly, in the resulting effective theory, the underlying 'classical geome- 
tries' will have a well defined Planck scale as defined by the expansion of the cosmological 
observers, only when the condition of boundedness of 9 is satisfied. 



IV. BIANCHI-I MODEL 

The loop quantization of Bianchi-I model is an interesting generalization of the isotropic 
loop quantum cosmology. Unlike the isotropic model, the Weyl curvature is not zero in this 
model and hence singularities are not only of Ricci curvature type but also of Weyl curvature 
type. As emphasized in our discussion earlier, this feature makes it necessary to go beyond 
the scope of energy density and its properties to understand singularity resolution. 

The spatial manifold in the Bianchi-I model can in general be non-compact with 
topology. Just as in the isotropic case, in order to introduce a symplectic structure, it is 
necessary to consider a fiducial cell V. With respect to the fiducial metric introduced on 
M^, this cell has volume Vo = hhh- The size of the fiducial cell and the fiducial lengths 
// are completely arbitrary freedoms which must be respected by the underlying theory 
at the classical as well at the effective and quantum level. As we will see the consistent 
implementation of these freedoms severely constrain the available quantization choices. If 
the spatial manifold is restricted to a compact topology such as a three-torus (T^), then the 
fiducial cell is no longer required as the integrations to define the symplectic structure are 
finite. However, we will show that even with this topology the choices which are restricted 
as above remain inviable. 
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The spacetime metric in the Bianchi-I model, for a particular choice of fiducial coordinates 
(x, y, z) is of the form: 

ds^ = -N^ dt^ + al dx^ + a\ dy^ + al dz^ (17) 

where iV is the lapse factor. The metric reduces to the isotropic one when the individual scale 
factors become equal. ^ Due to the underlying symmetries of the spacetime, the matrix valued 
connection and triad can be expressed as c/ and p/ respectively (where I = 1,2,3). 
These phase space variables satisfy 

{ci,pj} = 8nG^6jj (18) 

where 7 is the Barbero-Immirzi parameter. The metric components are related to the triads 
as 

bil = I2ha2a3, \p2\ = hkaias, {p^l = ^2^3 • (19) 

On the solutions of classical GR, the connection components are related to the time derivative 
of metric components as 

ci = -fliai. (20) 

It is useful to note that under the freedom of the change in the shape of the fiducial cell 
(^1,^2,^3) — ^ (^'15^2)^3)) (c/,P/) transform as 



1 



c 



I'lCi, Pi p[ = I'^l'^pi (21) 



and similarly for other components. 

Any quantization of this model should satisfy two consistency requirements: (i) If per- 
formed for a non-compact spatial topology, it must respect above freedom as well as the one 
in the choice of underlying fiducial metric. This corresponds to independence of physical 
predictions of invariant quantities such as expansion factor and curvature scalars from Vo 
and // and (ii) the quantum theory and its physical implications must reduce to those of the 
isotropic model in loop quantum cosmology when the anisotropic shear vanishes. The latter 
requirement ensures that the theory should have GR as its infra red limit at low spacetime 
curvatures. 

In the literature two different loop quantizations have been proposed which reduce to 
the consistent isotropic case in loop quantum cosmology. Both of these choices draw lessons 
from the 'improved quantization' of [4], albeit in a different way. The first choice relies on 
extending the prescription to regularize the field strength tensor in the anisotropic case "as 



is" in the isotropic case for triad variables [2J]. For the massless scalar field the resulting 



quantization can be shown to have certain similar features as the exactly solvable loo p q uan 



tum cosmology [10| and has been extensively studied in different contexts [25|, 126|, |27|, 128 
However, it turns out that this quantization suffers from lack of freedom for the choice of 
shape of the fiducial cell and one is therefore restricted to a compact topolo gy, such as a 
torus. The second choice for the regularization of the Hamiltonian constraint [22] (see also 



Appendix C of Ref. [25j where it was first discussed) successfully overcomes these limita- 



tions. At first sight it appears to be based on a prescription which is not as parallel to the 



^ There is still the possibility of changing the fiducial coordinates that will introduce a further freedom, in 
which case an isotropic geometry will not have the same factors a/. This fact introduces a subtlety in the 



passage from the Bianchi model to the isotropic one. For a detailed discussion of these issues see [22 1 
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improved quantization of isotropic LQC However, this viewpoint is incorrect. As we will 
show, this quantization in fact completely captures the 'improved dynamics' as expressed in 
/?, V variables. The first and second choice of regularization of the Hamiltonian constraint 
are sometimes referred in the literature as JI and fl' approach to reflect that the first choice 
mimics the JI or 'improved quantization' of isotropic LQC. Such a nomenclature is misleading 
for two reasons. The first is that in the isotropic limit both regularizations yield the JI 'im- 
proved dynamics'. Secondly, as we will show below, it is in fact the fl' approach which truly 
captures the regularization of Hamiltonian constraint as performed in isotropic LQC. How- 
ever, in order for comparison of our analysis with earlier works we shall use above notation. 
In principle it is possible to perform a general analysis incorporating other regularizations 
by introducing variables analogous to those in f|T5l) for the Bianchi-I spacetime. These will 
include those which do not lead to the unique consistent isotropic quantum dynamics in an 
appropriate limit. However, such models (which include the so called lattice refined models 
) suffer similar limitations as for the first choice of regularization (the JI scheme). These 



will not be considered in the following. 

We first discuss the quantization as proposed in Ref. [22i]. As we remarked earlier, the 
resulting theory respects the underlying freedom of the change in shape of the fiducial cell as 
well the fiducial metric. For the massless scalar field model, the quantum constraint has been 
shown to be non-singular for all states in the physical Hilbert space. As in the isotropic 
case, this model leads to a non-singular bounce of the scale factors when the spacetime 
curvature becomes close to the Planck value. For states which are semi-classical at late 
times, the underlying quantum dynamics can be expected to be very well approximated by 
the following effective Hamiltonian 

1 /sinfu'iCi) sin(7i;,c2) \ ^, 

^eff = - c ^ _/ ^ > iP2 + cyclic terms + T^^att 22 

where we have chosen lapse = 1 and 



P'l = ^a/-^' = ^a/-^' and /Xg = A./-^ . (23) 
V V V 

Above particular dependence on the triads arises from the regularization procedure employed 
in the quantization j22| . 

It is interesting to note that using eqs. ( !23l) one can rewrite the Hamiltonian constraint 
( I22I) in the following form 

^eflf = - o 9 X 9 ^ (sin(A/3i) sin (A/32) + cychc terms) + TYmatt 

where /?/ are conjugate variables to the volume V. In the limit A — 0, /^j ~ ^aijai. 
Note that the form of the constraint and the behavior of /3/ is analogous to the one in 
the isotropic case ([8]). Thus, as in the improved quantization of the isotropic LQC, this 
quantization implements a regularization of /5/ which is related to the expansion rate. 

An immediate consequence of (!22l) is that it leads to a universal bound on the matter 
energy density. The vanishing of the Hamiltonian constraint, Tiefr ~ implies 

P = o ^ 2A (sin(7iici) sin(7J2C2) + cyclic terms) . (24) 
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Since the terms in the parenthesis are all bounded functions, the maximum value of the 
energy density is given by pmax = 3/(87rGA^). Using Hamilton's equation for pj: 



Pi 



d 

inG-f—Hes 

OCi 



1 



cos(/i'iCi) (sin(/i2 C2) + sin(/i3 C3)) 



(25) 
(26) 



we obtain 

Pi ^ 
Pi 7-^ 

and similarly for the other terms. One can combine them and compute the expansion 9 of 
the congruences of cosmological observers 
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y_ 

V 



Pl ^P2_^P3 



Pi P2 P3 



27A 



[sm{fi[ ci + C2) + sin(/i'^ Ci + /i'a C2) + sin(/X2 C2 + /ig C3)] . 



(27) 



It is then obvious that the absolute value of the expansion | ^| is bounded, on the constrained 
surface as given by the effective dynamics. Just as in the isotropic model, it is expected that 
the spectrum of the corresponding operator in the quantum theory will be bounded. The 
maximum value of the expansion factor is given by 



1^1 



sup 



27A 



(28) 



This means that there is an invariant 'Planck scale' that binds this covariantly defined 
quantity and thus provides a well defined notion of a scale that appears as a true quantum 
gravity effect (note that it is inverse to the Planck length and thus would diverge in the 
'classical limit' ip\ ^ 0). 

Another consequence of the effective Hamiltonian (122|) is that the shear scalar 



a' := a'^V,, = - {{H, - H^f + {H^ - H,f + {H, - H,f) 



(29) 



also has a universal bound. This can be proved by finding directional Hubble rates Hi 
defined as the ratio hi/ai. A straightforward calculation leads to 



Hi — H2 



Ho — H-i 



H3 — Hi 



1 

7A 
1 

7A 
1 

7A 



sin(71^ci - /i2C2) + sin(773C3)(cos(772C2) - cos(7i'iCi)) 
sin(/i2C2 - /igCa) + sin(7iici)(cos(773C3) - cos(772C2)) 
sin(/i3C3 - Jl'iCi) + sin(772C2)(cos(7iici) - cos(7i3C3)) 



(30) 
(31) 
(32) 



Since each of {Hj — Hj) is bounded, it is obvious that cr^ is also bounded. 

It is important to note that for the effective Hamiltonian fl22|) . both 9 and are 
invariant quantities with respect to any freedom under the change of fiducial cell. In terms 
of these quantities it is quite transparent to reduce to the isotropic case, which corresponds 
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to dynamical trajectories for which Hi = H2 = H3, and therefore, o"^ = 0. Further, the 
upper bounds on both the entities are universal with values close to the Planck scale. 



Let us now consider the effective Hamiltonian which arises in the alternative quantization 
(the so called '/I' approach) ji^, 25, 26, 27, 28|. With lapse = it is given as 



H 



1 



eff 



87rG'72 



sin (/iiCi) sin {^2(^2) 
7^1/^2 



P1P2 + cyclic terms 



where 



A 



^2 



A 



and 



(33) 



(34) 



The vanishing of the above Hamiltonian constraint leads to the expression for the energy 
density: 



P = o ( — sin (piCi) sin {^2^2) + cyclic terms ) . 



(35) 



Even though the trigonometric terms in the parenthesis are bounded, p fails to have an 
upper bound because prefactors such as 1/^1^2/^3 can be unbounded. As an example, if 
the spatial manifold is a torus with a unit volume then ^Jpvpilvz = \/ai/a2 which diverges 
when ai ^ 00 or 02 ^ 0. Both of these limits are allowed possibilities in this model. Even 
when these prefactors are bounded, this bound depends on the phase space variables, as was 
the case for the isotropic case for a parametrization different from the improved one. 

To find the expression for the expansion factor 6, let us consider the equation of motion 
for pj. Using Hamilton's equation 



dpi Q 9 
dt oci 



(36) 



and dt = (l/l^)dt, we obtain 

dpi pi cosijLiCi] 



{P2 sin (/iaCs) + Pa sin (/igCg)) 



dt V-i 

and similarly for p2 and p^. The expression for 6 thus becomes 



(37) 



e 



1 



271/ 



Pl sin(/iiCi)(cOs(/i2C2) +C0S(;U3C3)) +P2 sin(/i2C2)(cos(;UiCi) +COs(/i3C3)) 



+ P3 sin (/i3C3)(cOS {piCi) + COS {P2C2)) 



(38) 



It is obvious that unlike the expansion factor in Eq. fl7r|) . this expression has no upper bound. 
Especially, the expansion factor diverges when either of the scale factors vanish (a behavior 
which is in common with the classical theory). On the other hand using (!27|) we note that 
for loop quantization of Bianchi-I spacetime as developed by Ashtekar and Wilson-Ewing 
' 2^ the expansion factor never diverges even when the volume vanishes. 
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The shear scalar (129|) shares the same fate as the expansion factor in this scheme. To see 
that we evaluate the directional Hubble rates from pi and find 



H^ — H'} = 



H2 — 



if 3 — Hi 



J93 sin (/X3C3)(cos {^2(^2) - cos (/iiCi)) + pi sin (yU^ci) cos {^1202) 



P2 sin {^12^2) cos (yUiCi^ 



Pi sin (/iiCi)(cos ifJ-^cs) - cos {^12^2)) + P2 sin (/i2C2) cos (/i3C3) 



- P3 sin (/i3C3) cos (/i2C2) 



P2 sin (/i2C2)(cOS (/iiCi) - COS (Ai3C3)) + P3 sin (/i3C3) COS (/iiCi; 



- Pi sin (/iiCi) COS (yU3C3) 



(39) 



(40) 



(41) 



Thus unlike the difference of Hubble rates in (!30|) . above expressions are unbounded implying 
that the shear scalar in loop quantization of Bianchi-I model 2J, |26|] may in general diverge 
for arbitrary matter content. 

This quantization prescription then shares the same properties that we analyzed in Sec. lIIII 
for the '/io quantization' (and lattice refining models). That is, even when there might be 
solutions to the effective equations of motion that 'bounce', there is no universal 'quantum 
gravity scale'. In the quantum theory this inadequacy manifests itself in the fact that the 
operators associated to 9 and will not be bounded. Just as we did for the isotropic 
case, we are led to conclude that this quantization, be it on M'^ or on the three-torus T^, is 
physically inviable. 

In summary, we have analyzed the effective theory one expects to get from two possible 
quantizations of the Bianchi-I system. We saw that the so-called Jl quantization does not 
yield a unique scale for the expansion and shear, whereas the /l' quantization does. Just as 
we did for the isotropic case, one can expect that the existence of bounded operators and 
therefore an invariant Planck scale [l3|, formulated in terms of the expansion and shear, can 
be taken as criteria for selecting the physically correct quantization. 



V. DISCUSSION AND OUTLOOK 

Let us summarize our results. We have put forward the proposal that, in order to address 
the issue of singularity resolution in any quantum theory of gravity, the natural geometrical 
quantities to consider for exploring the resulting quantum geometry are those quantities as- 
sociated to some observers. In the case of homogeneous cosmologies such natural quantities 
are given by the expansion and shear as measured by the cosmological comoving observers. 
The advantage of considering these observables is twofold. First, they posses a clear geomet- 
rical interpretation when there is an underlying classical spacetime, in terms of the behavior 
of the congruence of geodesies associated to cosmological observers. Second, since they are 
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well defined quantities, independent of any fiducial structures, they can also represent mean- 
ingful relational observables for the full quantum theory. They are the most natural with 
respect to which singularity resolution can be posed. Geodesies will be well defined and 
complete if these quantities are well behaved classically. Thus, one can expect that, if the 
quantum geometry is non-singular, the corresponding operators are to be bounded in the 
quantum theory. This also represent a step further from considering only the behavior of 



geodesies in an effective spacetime description as was done in [16 . 

With these geometrical considerations in mind, we have revisited the loop quantization 
of = FRW and Bianchi I cosmologies and focused on the issue of singularity resolution. 
A detailed study of these models lead us to argue that boundedness of these quantities in 
the quantum theory is a consistent requirement for singularity resolution, in our present 
restricted homogeneous context. We have seen that, if this condition is violated, as it 
happens in certain 'loop quantizations' for certain choices in the definition of the quantum 
theory, the resulting effective theory has the unphysical feature of not having a well defined 
'quantum gravity scale'. That is, the classical trajectories on phase space might represent 
universes that 'bounce', but such that the scale at which this happens is arbitrary, even at 
scales of ordinary phenomena, already ruled out by experiments. Another lesson we can 
learn from this result is that we can, in retrospect, see why the 'improved quantization' 
for the fiat isotropic model is the only physically viable choice. It is 'designed' to always 
bind the expansion, and therefore allow for geodesies -in the semiclassical scenario- to be 
complete. As we showed, this conclusion also holds true for the expansion rate and the 



shear scalar in Bianchi-I model for the Ji' quantization studied in Ref. [22|. Even though a 
quantization inequivalent to the one in Ref. 22j can be performed such that the theory leads 
to a correct isotropic limit j23 . 26 . 2?! . 28| . our analysis shows that it fails to be viable from 



the geometric perspective of singularity resolution. It is thus interesting that the uniqueness 
result for the isotropic model established in Ref. [l3| and in this paper (from a different 
perspective) seems to extend to the Bianchi-I spacetime. 

The criteria that we have put forward allow us to speak meaningfully about singularity 
resolution in the isotropic case, and extend it to the anisotropic scenario. In both these 
cases, out of the available quantizations, there is only one choice which leads to a physically 
viable resolution of singularities. Thus, the requirement of having a singularity free theory 
selects, in these cases, a unique consistent quantization. Still, even when one can be rea- 
sonably confident that this criteria is better motivated than the previous one available for 
the isotropic models |13|] , more work is needed to assure that it will continue to be the most 
reasonable condition for more general models such as inhomogeneous cosmologies. 

In this respect it is perhaps useful to end this contribution with a discussion regarding 
the issue of singularity resolution in quantum gravity. While it is generally expected that 
a quantum theory of gravity will 'cure' the singularities that are present in the general 
theory of relativity, it is not clear exactly how this is supposed to take place. Even in the 
restricted setting of homogeneous cosmologies, where only gravitational configurations with 
a large degree of symmetry are considered, this issue is not fully settled. For instance, if we 
take the viewpoint that in the classical theory (almost) all solutions possess a singularity, 
and if these classical solutions are to represent the phase space of the theory (as say, in 
geometric quantization) or be the dominant ones in the calculation of the path integral, 
then the resolution of the singularity by quantum effects seems somewhat puzzling. As is 
not so uncommon, reasonable expectations do not get realized when attempting to build 
detailed models for fundamental physics. Loop quantum cosmology has given an answer 



15 



that is somewhat unexpected. The way that the singularity is resolved is by changing 
drastically that set of 'classical configurations' that define the 'quantum phase space'. For, 
all solutions to the effective equations of motion are nonsingular, and those are precisely the 
trajectories that dominate the canonical path integral. A second point worth emphasizing 
is the way that singularities are resolved in LQC. It is not through the behavior of the 
wavefunction on the 'would-be singularity', nor through special boundary /initial conditions, 
but rather by analyzing the behavior of relevant Dirac observables. From this perspective, 
the main challenge is to identify those observables that are relevant in order to claim that 
singularities are generically resolved. As we have seen, in the homogeneous setting the 
choice is rather natural. Whether the criteria proposed here will still be meaningful when 
considering inhomogeneous cosmologies remains an open issue. 
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